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Abstract 

We have studied a quantum Hamiltonian that models an array of ultrasmall 
Josephson junctions with short range Josephson couplings, Ej, and charging 
energies, Ec, due to the small capacitance of the junctions. We derive a new 
effective quantum spherical model for the array Hamiltonian. As an applica- 
tion we start by approximating the capacitance matrix by its self-capacitive 
limit and in the presence of an external uniform background of charges, q x . 
In this limit we obtain the zero-temperature superconductor-insulator phase 
diagram, Ej rit (Ec,q x ), that improves upon previous theoretical results that 
used a mean field theory approximation. Next we obtain a closed-form ex- 
pression for the conductivity of a square array, and derive a universal scaling 
relation valid about the zero-temperature quantum critical point. In the lat- 
ter regime the energy scale is determined by temperature and we establish 
universal scaling forms for the frequency dependence of the conductivity. 



PACS numbers: 74.50. +r, 67.40. Db 
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I. INTRODUCTION 



Quantum phase transitions have attracted a significant amount of interest in recent years. 
There are different physical systems where Quantum phase transitions can be studied, like in 
magnetic systemsQ, the quantum Hall effectll, superconducting films! and Josephson junction 
arrays (JJA)i. There have been significant advances in lithographic techniques that allow 
the fabrication of arrays of ultra-small superconducting islands, with charging energy, Eq, 
that can dominate the Josephson coupling energy, Ej, making quantum fluctuation effects 
of paramount importance!. Because the junction parameters can be controlled accurately 
JJA offer a unique system where one can test the nature of quantum phase transitions and 
critical phenomena, in particular the superconductor-insulator (SI) phase transition induced 
by quantum fluctuations. Therefore JJA can be a prototype system that can display well 
controlled quantum SI phase transitions. 

In the JJA problem, the two relevant low temperature energy scales are Ej, that permits 
tunneling of Cooper pairs between islands, and Eq that tends to localize the charge carriers 
in the islands. When Ej is much larger then Eq the phases in the islands are well defined. 
In this regime only the phases determine the properties of the JJA. In the opposite limit, i.e. 
Eq » Ej the superconducting phase is perturbed by strong zero point quantum fluctua- 
tions due to the Coulomb blockade that localizes the charge carriers to the islands. Several 
theoretical HI and experimental studies! have considered the competition between the Ej 
dominated phase and the E c dominated charging energy regions in periodic JJA. It has been 
established that for sufficiently large charging energy the quantum phase fluctuations lead 
a complete suppression of long-range phase coherence even at zero temperature. This type 
of quantum phase transition has attracted significant interest in recent years (for a review 
see Ref. |H] and also |XT|) . Generally, these transitions take place at zero temperature where 
crossing the phase boundary is accompanied by a change in the ground state of the system. 
This transition is induced at zero temperature by changing some external parameter in the 
Hamiltonian of the system, for example the charging energy in the Josephson-j unction ar- 
rays (given by the quantum parameter a = Ec/Ej). Let the critical temperature be T c and 
denote the distance from criticality as 5t = 1 — T/T c . The critical behavior is asymptotically 
close to the critical point which is entirely determined by classical physics. This is because 
the characteristic frequency associated with the critical fluctuations, uj c ~ 1/£ T) vanishes at 
criticality, and the characteristic correlation time for relaxation towards equilibrium £ T ~ £ z 
diverges at the critical point (here £ ~ \t\~ u is the diverging characteristic correlation length 
in the system.) A quantum system behaves classically if the temperature is larger than all 
frequencies of interest, i.e. as long as hu c « k B T, which is always the case when the system 
is sufficiently close to the critical point with a non zero T c . In the critical region various 
physical quantities can have scaling forms (i.e. given in terms of homogeneous functions) 
as a consequence of the diverging £ and £ T . Since the experimentally accessible systems are 
necessarily at non-zero temperature, one needs to understand how the quantum behavior 
is modified at finite temperature. The important observation is that at finite temperature 
the partition function has a finite temporal dimension, since the effective classical system 
extends to an extra dimension of size /3h. 

At T = 0, the 2D JJA problem goes into a 3-dimensional classical XY model with the 
critical behavior of the 3D XY model. Moreover, there should be a clear signature in the 
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nature of the correlation functions (like for example in linear response) when crossing over 
from the 2D XY model at high temperatures to the 3D XY model as T — > 0. Specifically, 
the response of the quantum critical point to a small temperature perturbation will be 
determined by the location of the quantum-to-classical crossover. An especially interesting 
quantum critical regime appears when the relaxational frequency satisfies Tluj ksT. In 
this regime the only relevant scale is given by the temperature and the dynamical scaling 
functions will only depend on the ratio hu/kBT. As the ratio hu/kBT is varied in an 
experiment we expect to see a crossover between temperature- and frequency-dominated 
scaling regimes. For the finite temperature experiments one can scale the data by using the 
knowledge of the frequency dependent scaling functions to test the existence or probing of 
a quantum critical point. The experimental identification of a quantum phase transition in 
JJA will rely upon finding the scaling behavior with the relevant parameters as temperature, 
frequency and the wavelength dependence of various observables. The specific signatures are 
given by universal values of certain dimensionless amplitudes as, for example, the resistivities 
at the critical point in 2D JJA systems. Motivated by these issues, it is the goal of the present 
paper to study the superconductor-insulator transition in 2D JJA in the quantum critical 
regime, by focusing on the scaling properties of the directly measurable electromagnetic 
conductivities in JJA systems. 

The outline of the rest of the paper is the following: In Section II we define the model 
Hamiltonian, followed by a novel path integral formulation of the problem. In section III we 
present as a test of our quantum spherical model approach for the self-capacitance model. 
There we show that the spherical technique gives a better answer than the one obtained 
from mean field theory. In Section IV we present our calculation of the scaling properties 
of the conductivity in the quantum critical regime. Finally, in Section V we briefly discuss 
our results. 

II. THE MODEL 

A. Quantum phase Hamiltonian 

A Josephson junction array (JJA) can be modeled by a periodic lattice of supercon- 
ducting islands separated by insulating barriers. Each island becomes independently super- 
conducting about the bulk transition temperature T c0 and it is characterized by an order 
parameter ip( r i) = l^o^i)! 6 *^ 1 "^! where is a two-dimensional vector denoting the position 
of each island. The magnitude of the order parameter, |-?/>o( r i)|> is non-fluctuating when 
the temperature is lowered further and the onset of long range phase coherence due to the 
tunneling of Cooper pairs between the islands is responsible for the zero resistance drop in 
the arrays. The phase coherence onset temperature in JJAs can be significantly reduced 
by making the junctions ultrasmall. When the junction's capacitance is small the charg- 
ing energy (i.e. the energy necessary to transfer Cooper pair charges between the islands) 
increases to the point where no pairs can tunnel anymore, completely quenching the Joseph- 
son current. The competition between the Josephson tunneling and the charging (Coulomb) 
energy, without dissipation, can be modeled by the Hamiltonian 

H= H c + Hj, 
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He— — — J^[C l } rr 'Q r Q r >, 

A r 

Hj= £ J(|ri-r 2 |)[l-cos(«/» ri -0 r2 )]. (1) 

( r l: r 2> 

Here, Qr = (2e/i)d/d(p r is the charge operator while r represents the superconducting phase 
operator of the grain at the site r; J(|i"i — r 2 |) is the site-dependent Josephson coupling and 
[C _1 ] rr ' is the inverse capacitance matrix. 



B. Euclidean action and path integral representation 

It is useful to derive a field-theoretic representation of the partition function for Eq (1) 
to study the quantum nature of the superconductor-insulating (SI) phase transition. A 
convenient procedure is to introduce a path-integral representation in a basis diagonal in 
4>j. In this representation the partition function is expressed as 



e -Sl0J. ( 2 ) 




Here the functional integral is evaluated over the phases restricted to the compact interval 
[0, 27r] and with an effective action (% = 1) 

SM = S c [0]+S,M, 

r,r' 

SM = J^dr £ J(| ri -r 2 |){l-cos[<MT)-0r 2 (T)]}. (3) 

" ' < r l. r 2) 

Because the values of the phases which differ by 2tc are equivalent, the path integral 
can be written in terms of non-compact phase variables 6 l r (r), defined on the unrestricted 
interval [— oo,+oo], and by a set of winding numbers {n r }, which are integers running 
from — oo to +oo. Consequently, the path-integral representation in Eq.(||) also includes a 
summation over winding numbers^: r (r) = 8 r (0) + 2Trin r r/P + 9 r (r) which reflects the 
discreteness of the charge, so that the integration measure reads 

^ = e n / ^(°) / vd M- ( 4 ) 

fcyrJo Mo) 

This concludes the formulation of the model. In the following section we will transform 
this representation to one where the order parameter field is expressed in a novel quantum 
spherical model self-consistent scheme. 



III. QUANTUM SPHERICAL MODEL APPROACH 

To study the JJA model it appears at first natural to use a description in terms of 
an effective Ginsburg-Landau functional derived from the microscopic model of Eq. 
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Several studies of JJA have followed this route, also known as the coarse grained approach 
first developed by Doniach0. The key point of this method is to introduce a complex field 
order parameter i/? r (or equivalently a two component real field) whose expectation value is 
proportional to the S r (0) vector defined by 

S r (</>) = [^(0),^(0)] = [cos(0 r ),sin(0 r )]. (5) 

The non-zero thermal average (S x (<p)) = (cos (pi) describes the "phase-locking" or long-range 
phase ordering in the model. The system is governed by the Ginzburg-Landau functional 
as long as the order parameter is small and the decoupling of the Josephson energies is 
valid. This is a serious shortcoming of the coarse grained approach since the method is 
restricted to the region of parameters in the vicinity of the critical point and does not offer 
a self-consistent description of the full problem. 

There are no nontrivial exact solutions of the model, except for some molecular field- 
type approaches. It is therefore reasonable to develop new approximate mappings of the 
model which may admit exact solutions, both in the ordered in disordered phases. ^From 
the definition of the pseudo-spin variables S r the following rigid constraint holds for each 
one, 

|S r (0)| 2 = [^(0)] 2 + [^(0)] 2 = cos 2 (0 r ) + sin 2 (0 r ) = 1. (6) 
The relation in Eq. also implies that a weaker condition also applies, namely: 

E|Sr(0)| 2 = iV. (7) 
r 

The main idea of our approach is to attempt to generate an effective partition function 
from the original one with cosine interaction, which incorporates the constrained nature 
of the original variables. This leads us to the formulation of the problem in terms of the 
spherical mo del&H by introducing the appropriate constrained order parameter field. The 
name of the model comes from the observation that in Eq.(^) the allowed states of the 
spherical model are all points on the surface of a hypersphere of radius N. Although the 
spherical model was originally introduced as an approximation to the classical Ising model, 
the "sphericalization" technique may be readily applied to a variety of other problems (for a 
quantum generalization of the spherical model, see Ref. |20|); the one key necessary ingredient 
is the existence of an inherent constraint on the Hamiltonian variables as in Eq.(gj). The 
model defined by Eq.(^) is in fact a hybrid of the spherical and the two component (M = 2) 
vector models. The global constraint of Eq. (^) may be introduced on the continuous order 
parameter field by using the functional analogue of the Dirac-delta function 



Hj2\^r(r)\ 2 -N 



IP [Re^ r (r) - S?(0(r)] 5 [Imt/> r (r) - 5»(0(r))] , (8) 

r 

where the i/'r(' r ) are the complex c-number fields, which satisfy the quantum periodic bound- 
ary condition if) r (/3) = il> r (0), and taken as continuous variables, i.e., — oo < V , r( r ) < +oo, 
but constrained (on the average) to have unit length. The partition function of Eq. (Q) then 
reads 
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JV e' 



X 



T 



x5[Im^ r (r)-^(0(r))]. 



(9) 



One may wonder whether the introduction of the spherical condition on the order parameter 
fields in Eq. (|7|) is essential since the global constraint is automatically fulfilled by virtue 
of the exact relation (^). Indeed, in a rigorous treatment of the partition function of Eq. 
(0), the constraint will be spurious; in any approximate treatment of the functional integral 
for Z, however, (as for example in a coarse-graining approach) the introduction of the 
unrestricted order parameter field generally will lead to a violation of Eq.(^). In such a case 
the introduction of the spherical constraint on the order parameter field in the partition 
function simply reflects the constrained nature of the original pseudo-spin variables. 

Integrating over the phase variables leaves the statistical sum entirely written in terms 
of the constrained continuous order parameter field t/) T (r), so that the partition function 
becomes 



^El^(r) 



N e 



-Sefl#] 



where the effective order parameter action reads 

5eff[^]=5j[^] + $ C [V>], 

with 



(10) 



(11) 



(12) 



<ri,r 2 > 



and 



-*c#] 



n 



2ni 



cxp 



(j2^dTfl r (T)-i> r (T)-W[ri 



(13) 



where the generating functional of the cumulant multipoint phase correlators of the "non- 
interacting" system (i.e., involving only charging energy terms) is 



W[fi] = \nj n^0r exp J* dr^{r) ■ S r (0(r)) + <S, 




(14) 



with the variables \x acting as the source fields. It can be seen that from Eq.flTjj) that 
e -$a[ij>] j s j us t the functional Fourier transform of e~ w ^. The standard way to proceed is 
to calculate $c[' 1 /'] using the saddle point method and a subsequent loop expansion in terms 
of the powers of the order parameter using Eq. (|T4]) . The structure of this expansion is 
briefly described in the Appendix. 

It is convenient to employ the functional Fourier representation of the 5-functional to 
resolve the spherical constraint in Eq. (|i0|) : 
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Accordingly, we can write 



+ioo 



2m 



x exp 



dri{r)g(r) 



(15) 



""J 

x exp 



+ioo 



VX[r) 
2m 



x 



[dr\{r) (j:\^r(r)\ 2 - N 



(16) 



We may factor the trace over each iI> t {t) at the expense of introducing a new integral over 
one component field A(r). Expanding part of the effective action $c[V'] to second order in 
i/> r ( r ) we get from Eq.(|iO|) 



V QSA 



vx 

2~ni 



(17) 



where 



Sqsa[^,A]= jfVdr'j 2 [(J(ki-r 2 | + A(r)<J ri ^)(J(r-r') 

"" L <ri,r 2 ) 

+ ro+^r; r 2 r')] <(t)^ 2 (t') - iVA(r)<f(r - r')} , 



(18) 



is the effective action of the quantum spherical approximation (QSA). Subsequent improve- 
ments over the spherical model approach may be obtained by considering the saddle point 
corrections to the Lagrangian parameter A. This can be conveniently done in terms of a loop 
expansion as described in the Appendix. Furthermore, ro 2 + (rir; ^t') = [W^ 1 ] h ( r i r ; r 2T / ) 
is the two-point phase vertex correlator and 



W, 



02 '{tit;t 2 t 



with 



Z ° {n r } r 



^ = En 

{n r } r 



2tt 



d6(Q) 



6»(0)+27rn r 
9(0) 



(19) 



2jt 



d0(O) 



6»(0)+27rn r 
0(0) 



m.Me- 5c[ ' 



(20) 



is the partition function for the "non-interacting" system. Eq. fllBD incorporates the 
quadratic term proportional to A(r) and the integration in Eq.([l7|) takes place over all 
configurations of the order parameter field. The model is now unconstrained and quadratic, 
so all quantities can be readily computed. In the thermodynamic limit, iV —>■ oo, we can 
calculate the functional integral in Eq. ([T7J) by the steepest-descent method. To proceed, we 
introduce propagators associated with the order parameter field defined by 

G(r 1 r;r 2 r') = «(r)^ r2 (r')) QSA . (21) 

The condition that the integrand in Eq.(|l7|) has a saddle point A(r) = Aq is that 
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l = ^£G(rT;rT + + ), (22) 

r 

which becomes an implicit equation for X . The QSA ensemble average is now defined by 

/[n r ^ r ^]---e-^^A ] 
QSA " J[U r ViP r V^ r }e- s ^M] 



A Fourier transform of Eq. flT8| ) in momentum and frequency space enables one to write the 
spherical constraint (|2"2] ) explicitly as 

1 = _L i f 24) 

^ k £ A- J(k) + [^ 02 (^)]- 1 ' 1 ; 

with J(k) the Fourier transform of the Josephson interactions J(|ri — r2|), W^i^e) the fre- 
quency transformed phase-phase correlator of Eq. ([19]) and ug = 2n£/ j3 (£ — 0, ±1, ±2, . . .) 
the (Bose) Matsubara frequencies. As usual in a spherical model analysis the critical behav- 
ior is determined by the denominator of the summand in the spherical constraint equation 
of Eq. (p2|) . Specifically, when [1/G(k = 0,ue = 0)] = 0, where G(k, ujg) is the Fourier 
transformed order parameter correlation of Eq. (|2T|), the system displays a critical point at 

A - J(k = 0) + [W 02 {ut = 0)]- 1 = 0, (25) 



provided the momentum and frequency summations in the constraint in Eq. (|22[) converge. 
The system does not show a phase transition if the frequency-momentum sum in Eq.(f22"D 
diverges for N — > oo at the point defined by the criticality condition (p5|) . The universal 
critical properties only depend on the low-frequency behavior of [Wo2(<^)] _1 and the long- 
wavelength properties of the interaction J(k). Formally, the T = critical behavior will be 
identical to that of interacting quantum rotors with a kinetic energy proportional to k 2 in 
d = 3 dimensions, i.e. the transition will be in the universality class of the tree dimensional 
(3D) XY model. 



A. Zero temperature ground capacitance model. 

Offset charges are an important ingredient in the experimental array samples made of 
ultrasmall junctions. In this subsection we reconsider this problem as a test for the usefulness 
of the quantum spherical model approach. The offset charges, or an external gate voltage 
applied between the array and the substrate, behave like a chemical potential for injection 
of Cooper pairs into the array. Several authors have shown that static background charges 
can have a pronounced effect on the SI transition at zero temperature0'0. Including the 
offset charges, q x , in the charging energy of Eq. @ gives 

5 *]^r-DC-v(f-.)(^-.). (26) 

Furthermore, we simplify the model to include only the background capacitance (or self- 
charging) model. In this case [C _1 ] rir / = 5 r y/C and E c = |e 2 [C _1 ] rjr = e 2 /2C*o. Of course 
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the approximation that Cq » G\ y where G\ is the mutual capacitance between the islands, 
is a first step in the analysis that leads to physical insights into the general problem. Finally, 
we assume a square array characterized by the nearest-neighbor Josephson coupling Ej with 
J(k) = Ej[cos(k x a) + cos(A; y a)], 

with a the lattice constant. We obtain the corresponding density of states 



K 



\e\), 



(27) 



where K.(x) is the complete elliptic integral of the first kind0. The phase-phase correlation 
function becomesS 



W 02 (uji) 



SE C ^ ex P [-4pEc{q ~ Qx) 
~Zo~ Y {±EcY ~ [8E c (q - q x ) - iu t 



2 ■ 



(28) 



where Z = J2 q ex P [~ ^/3Ec(q — q x ) 2 ], and the summation is performed over all integer- 
valued charge states q — 0, ±1, ±2, . . ., which makes the function W 02 (^i) periodic. At low 
temperatures the sum over q in Eq.(|28D is dominated by the charge q which makes the 
exponent in the numerator of Eq.(p8|) smallest. For T = this value is q = and the 
equation for the critical line (E/Ec) vs. q x is obtained from the implicit equation 



1 



dE- 



p{E) 



\ + 2E c -E- 8E C (q x + ^) 



2 • 



(29) 



valid for — | < q x < | (other values of q x are included by a periodic extension with the 
period equal to one). The criticality condition (cf. Eq. (|25D ) reads 



By formally setting Ao 
solution: 



Ao = 2Ej - 2E C + 8E c q 2 x . (30) 
in Eq.(RDJ) we obtain the coarse-grained mean-field parabolic 



Ej_ 

Ec 



l-4^ 2 . 



(31) 



By substituting the value of Ao from Eq. (|30D into (|^) within the QSA method, after per- 
forming the summation over Matsubara frequencies, one obtains the T — > limit result 



1=P dEp(E)^ 



E, 



c 



2(A + 2E C - E) 



Ec 
: Ej 



sign | 4 ?x^ 



^2(A + 2E C - 


■E)E C \ 


) 


yj2(X + 2E C 


-E)E C 



E, 



(32) 



where V denotes the principal value of the integral. After Eq. 
phase boundary of the insulating Mott lobe at zero temperature 



we finally obtain the 
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1 - 2/ 2 /4 



(33) 



In Fig.(|]) we plot the resulting phase boundary in the (Ej/E c )-vs-q x plane at T = 0. 
We recognize the periodic lobes in q x of the insulating phase separated by regions of phase 
coherent superconducting state. For q x = 1/2 the superconducting state extends down 
to arbitrarily small values of Ej/Eq- We can compare the phase diagram resulting from 
the QSA to the one obtained from mean-field theory via the coarse grained approach (cf. 
Eq.(pTD). Whereas QSA gives the value Ej(q x = 0)/Ec ~ 1.653, mean-field theory gives 
Ej(q x = 0)/Ec = 1, which underestimates the critical value of Ej. It is also instructive to 
compare the QSA result for Ej(q x = 0)/E c to the recent third-order perturbation expansion 
in Ej/EM 



1 
2 



lEj_ 
±E C 



Ej_ 
128 \E C 



11 (E 



1025 \E C 



+ 



El 

e c 



(34) 



which gives Ej(q x = 0)/Ec ~ 1.590, giving a 3.8% of the difference from the critical value 
of Ej/Eq at q x = obtained in the QSA. 



IV. CONDUCTIVITY SCALING 

The conductivity is an experimentally measurable quantity in JJA. There are other 
studies of the superconductor-Mott-insulator and its universal conductance in JJA. Cha et 
al (see, Ref. [H]) carried out a 1/N expansion and Monte Carlo analysis. Van Otterlo et. al. 
used the coarse-graining approachll and Fazio and Zappela did an e-expansionfU. Here, we 
are interested in another aspect of the SI transition, namely, the scaling of the conductivity 
in the vicinity of the quantum critical point. 

To evaluate the conductivity, we need to add an external field in terms of a minimally 
coupled order parameter to the vector potential A(r, r). The Josephson coupling term in 
Eq. fl3|) then becomes 

J(|r, - r,|) -> J(|r, - r,|) exp £ A • dl) . (35) 

The phase shift on each junction is determined by the vector potential A of the magnetic 
applied field and in a typical experimental situation it can be entirely ascribed to the external 
field. 

The imaginary-time frequency dependent linear-response conductivity is given by 

<vK,q) = ^-j<? r [ ^/ r ;ii (0 ,or T+,q - ^ 

For vanishing magnetic field and offset charge the longitudinal component of ct(uj u ) = 
a xx {uj u , q = 0) in the QSA is 
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a(u v ) 



dp 2 F, 2 rir/a rib rlk 1 



%L0 V J-n/a (2n/a) 2 (5 



(37) 



where G _1 (k, u)e) = [A — J(k) + 2E C + u 2 /8E c ], and A is determined self-consistently from 
the constraint equation ([?[). To proceed it is convenient to obtain the generalized density of 
states for the 2-D square lattice as 



P2(E) = / 

J —IT 



7r / a dk x dk 
la (2n/a) 2 



7T 2 E, 



v sm 2 (ak x )5(E-J(k)) 

2 



E 



H- 



E 2 
AE> 




Q (AEj — E 



(38) 



where E(x) is the elliptic integral of the second kind. In terms of Eq. 
of Eq. (|37j) becomes 



the conductivity 



a{u v ) 



Ae 2 E 2 
Tiuj u 



/+oo 1 
dE- J2 P2(E)G(E, u t ) [G(E, w t ) - G(E, lu u + u t )] . 
-°° P Lot 



(39) 



Here, G~ 1 (E, u t ) = [5\ - E + 2Ej + uj/8E c ], where 5 > 



A- A 



crit- 



The zero-temperature 



critical boundary of the phase coherent state is signaled by \(Ej/Ec, T) reaching the value 
X cr it(Ej rit / Ec , T = 0). The parameter 5\ plays the important role of energy scale (which 
vanishes in the superconducting phase). We now consider the T = phase transition 
between the long-range ordered coherent phase state and the insulating phase by varying the 
coupling constant Ej through the critical value EJ l1 , where there is a diverging correlation 
length £ ~ \Ej — -Ej nt |. At finite temperatures, the deviations from T = behavior are 
distinguished by the scale set by the thermal coherence length £^ ~ T~ l l z (where z is the 
dynamic critical exponent). The quantum-critical region is defined by the inequality £r < £. 
In this case the system feels the finite value of the temperature before becoming sensitive 
to the deviations of Ej from E c j lt . In this regime the dynamic conductivity turns out to be 
remarkably universal. To proceed, we need to determine the 5\ dependence of temperature 
and Ej/Ec- The spherical constraint of Eq. ([?[) takes the form 



iv 



dE- 



P{E) 



5 X + 2Ej -E-SE C (4^-) 



2 ■ 



(40) 



\8E C J 



The near-critical properties of the spherical model are essentially determined by the struc- 
ture of the spectrum of the interaction matrix J(|i"i — r 2 |) in the neighborhood of its upper 
boundary, specifically by the behavior of the density of states associated with J(|ri — r 2 |). 
We therefore expand the square lattice density of states (^7]) about the upper limit of the 
spectrum of J(k): 



1 



EjTC 2 



\ 



(E-2Ej 
4Ej 2 



1 



+ 



E 



2Ejtt 8E/ir 



+ 0(E 2 ). 



(41) 



The subsequent summation frequency yields 
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2vr E 



In 



J 



sinh 




'2£ c (4£./ + 5 A ) 


sinh 


^2E C {8 X )_ 





(42) 



This constraint equation can be explicitly solved for the dependence of 5\ on T, Ej, and 
giving 

(k B T) 2 



2E r 



arsinh 



e -WEj sinh pJ 8 E c E, 



and near the T = quantum critical point 



(k B Tf 



2E, 



< arsinh 



c 



\ exp (-7t{3(Ej - Ef *))] } 



(43) 



(44) 



where i?} rit = 2Ec/tc . We extract the conductivity, as a function of the real frequencies 
lj in the form of real and imaginary parts of <r{uS) = cr Bins (uj u ) + a reg (u u ), after performing 
the sum over Matsubara frequencies in Eq.(^) followed by analytic continuation to real 
frequencies iue — > u + i0 + getting 



Imff S mg(w) 



Imcr reg (a;) 



and 
Re<x 
Recx 

x5 



4e 2 E 2 
' h 

Ae 2 E 2 



+oo i 8/3E c cosech 2 (f ^8E C (5 X - E + 2E. 

-oo dEp2 ^Z 5 A - £ + 2£j 



d£p 2 (£) 



1 



1 



to 



coth (§^8£ c (5 A -£ + 2£j 



4\/ 8E c (5 x -E + 2Ej)3 



SE C 



4(5 A -E + 2Ej) 



(45) 



sing 



UJ 



Ae 2 E 2 it 

' n 

2 ■ 4e 2 E> 
ft 



dEp 2 {E)5{uj 
1 



dEp 2 (E) 



8/3£ c cosech 2 (£yj8E c (5 x - E + 2Ej 
5a - £ + 2£j 

'0 



4 (<5 A - £ + 2E,j 



coth Z-J8E c (6x-E + 2Ej) 



UJ 

8E~c 



- 4(5 A - £ + 2Ej 



(46) 



The real part of the conductivity contains two contributions: the first, Recr s i ag (a;), is singular 
since it is proportional to 6(u) and the second finite-frequency is regular, Re<7 reg (u;), which 
arises from the electromagnetic field induced transitions to excited states. The singular part 
in turn is due to the free acceleration of quasiparticles. This is so since the JJA model 
considered here contains no dissipation mechanism which would e.g. arise from a coupling 
of the phase degrees of freedom to the normal electrons (Ohmic damping). At T = the 
singular part vanishes while the regular one can be evaluated explicitly by performing an 
energy integration in Eq. fl46|) with the result 



Re<7 reg (a;, S) 



(2e) 2 1 
h 25 



UJ, 



UJ 



T 



6 



where J-(x) = E [\fl 



X' 



x 



1 + 



UJ 



e 



UJ 
UJ C 



1 Jl+A/5 



UJ 
UJr. 



, (47) 



x- 



UJ 



uj/uj c , with u c = y/32E c 6x and 6 



5\/Ej. Here u c denotes the thrashed frequency above which the particle-hole excitations 
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can be created and the real part of the conductivity is finite, while Rea ieg (uj) vanishes for 
uj < u c indicating a Mott-insulating phase. Letting the gap parameter 5 go to zero, while 
keeping the ratio uj/uj c finite, and using the result 



lim — 



x = — 



(48) 



we obtain from the general QSA result (E7|) (valid for arbitrary distance 5\ away from the 
critical point) the long-wave length limit near-critical form of the conductivity: 



lim Re<7 reg (ci;, 5) 



(2e) 



7T 



<5^0 



h 




UJ 
0J r 



- 1 



(49) 



A plot of this result for different values of 5 is shown in Fig. (|2|), where we clearly see an 
asymptotic gap as we approach the critical point at zero temperature. This analytical result 
was previously derived in studies that relied on the coarse-grained and loop-expansion ap- 
proaches. From Eq.(|49|) at the transition, where the thrashed frequent uj c vanishes, a finite 
dc conductance cr* = (tt/8) ■ 4e 2 /h emerges which is the universal zero temperature conduc- 
tivity earlier obtained by Cha et al. In practice, however, all experiments are performed at 
low but nonzero temperature. 

We will now present the scaling analysis satisfied by <j(uj) in the vicinity of the quan- 
tum phase transition Ej = Ej rit . The temperature is taken to be nonzero but must obey 
k B T « Ej. A nonzero T implies the absence of long-range phase coherence and the scaling 
properties will depend upon a variable which measures the distance of the superconduct- 
ing ground state from criticality. The behavior of the conductivity in this regime can be 
understood in terms of the scaling forms 

R , x (2e) 2 ,( uj (k B T) l/vz \ 
Recx(c<;) = — ; — g 1 1 



Imer(a;) 



h 

(2e) 2 
h 



-9 



k B T' SEj j 
uj (k B T) l l vz * 
k B T' 6Ej 



(50) 



where g'(X,Y) = g' aing (X, Y) + g' Ieg (X, Y) and g"{X,Y) = g" 9ing (X } Y) + g" rcg (X, Y) are 
highly non-trivial but universal two-parameter functions. ^From equations ([45]) and ([46]) 
and the explicit solution (03) of the spherical constraint equation we derive 



*W*, Y) = \5{X)~\Y) du (u - 1) s(y, 

E^Y) 1 



u 



d'reAX, Y) — — 



1 - 



X 5 



cotli [ — ]0 



X< 



~\Y) 



(51) 



where 



E(Y) = 4 In Q v / e -W^ + 4 + ^e- 7r/y ^ 
Z(Y,u) 



+ 4 + \e-*l y ) U - (|Ve- 2 ^ + 4 + \e-l Y 



2 ! 



(52) 
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and 



Y) = j™du(u- 1 -) u) 



1 X M Q(Y,u) 



with 



Q(Y u) = V - - V - L— (54) 



We are in the quantum-critical regime by setting SEj to zero (Y = oo), while keeping the 
temperature small but finite (X ^ 0), with z = 2, and zv — 1, with the critical fluctuations 
quenched in a universal way by temperature. Here, k^T appears to be the dominant energy 
which determines the physics of the problem. At large frequencies, hoo » ksT, the finite 
temperature effects do not become manifest and the system displays the behavior of the 
zero-temperature critical point (5Ej = 0,T = 0). However, at small enough frequencies, 
Tiu) ~ k B T, the finite temperature effects of a become apparent which introduces a new 
energy scale into the problem. While fluctuations with frequencies larger then kgT '/h are 
unaffected those with hu < k^T will behave classically. We show the regular parts of the 
nontrivial scaling functions g'(X,Y) and g"{X,Y) in Figs. (W) and (HI). 



V. SUMMARY 



In this paper we have developed a microscopic analysis of a coupled array of ultrasmall 
Josephson junctions in the presence of charging energy effects. Our analysis is based on the 
quantum spherical model approach and the path-integral formulation of quantum mechan- 
ics explicitly tailored for the microscopic JJA Hamiltonian. The effective action formalism 
allows for an explicit implementation of the Coulomb and offset voltage effects into our 
consideration. Using this formalism we have considered the zero-temperature phase transi- 
tion to try to understand the finite temperature behavior of the system in terms of scaling 
functions for the electromagnetic response of the array. This is important since the exper- 
imental analysis of quantum phase transitions relies on the scaling behavior of observables 
(eg. the conductivity) for relevant physical parameters (like temperature and frequency). A 
relevant question which naturally arises is: To what extent does this phenomena depend on 
the model used to describe the JJA? Future theoretical studies on quantum critical behavior 
should consider disorder effects (random offset charges), coupling to quasiparticles (dissipa- 
tion), magnetic field frustration and finite size effects, just to call a few relevant issues in 
this problem. The formalism presented in this paper can be used to look at these specific 
situations. 
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APPENDIX A: LOOP EXPANSION FOR THE EFFECTIVE ACTION 



The generating functional of Eq.([T3|) for the connected cumulant functions 
) is given by the functional expansion 



00 1 r 

= V — / dxi . . . dx m W 0m (xx, x m )n(xi) . . . n{x m ), 
where we have for convenience introduced the short-hand notation 



(Al) 



to obtain 



J dx... = dr m X] X] • • • ' 



Wom(xi, ...,x m ) = (S%[<Kn)] • • • S a r zW m )]) c r- 



(A2) 



(A3) 



|cum 



For a given set of operators A,B and C the cumulant averages are defined as: (AB)[ 
(AB) -(A) (B) ; (ABC) C = (ABC) - (A) (BC) - (B) (AC) - (C) (AB) +2(A) (B) (C) , 
. . . etc... where 



(•••)( 



I[UrV(j> r } e-scl 



(A4) 



which is closely related to the multi-point vertex function T 0m (xi, • • • , x m ). The generating 
functional reads: 

00 1 f 

r [^] = Y] — r / dxi . . . dXmTomiXi, Jc m )V>(zi) . . . if>(x m ). (A5) 

The relation between W[/j] and T[i/j] is then given by the Legendre transform 

F[if>] + W\fjt] = J dxfji(x)il>(x). (A6) 

We may now evaluate the action $c[V'] defined in Eq. (|l"3"D , by introducing in the exponential 
of Eq. (|13"D a formal parameter ( to order the perturbation expansion (( will be set to unity 
at the end of calculation). We obtain 



^ = / [n^] ex P [-C fe jfWO") • tf r (r) - 



(A7) 



The steepest descent procedure is the standard method to obtain the expansion for $c[i/?]. 
The first step is to determine the saddle point /x of the integrand in Eq. ( |A7|) followed by a 
systematic expansion in terms of fluctuations about the saddle point. The procedure is by 



now standard (see, e.g. Ref. El) and we quote the final result: 
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cxp 



X 1 + 



J dxi . . . dx^W^{xi, . . . , X4, i/))D(xi,x 2 ; ij))D(x 3 , x 4 ; 
+ J dx l ...dy 1 W {3) (x 1 ,x 2 ,x 3 ;ij>)W (3) (y 1 ,y 2 ,y 3 ;ip) (^D(x 1 ,x 2 ; xj)) x 
x D(y u y 2 ; ip)D(x 3 , y 3 ; + ^p{x u y x \ i/>)D(x 2 , y 2 ; if>)D(x 3 , y 3 ; iff) 



-2^ 



(A8) 



where 



<W[/Lt] 



5(jl{x) 



+ ±- f dy 1 dy 2 W®{x,y 1 ,y 2] il>)D(y 1 ,y 2] il>) + 0(C 2 ), 

MoO^) = ^0*0 + ^J d y d yi d y2W {3) (x,y 1 ,y 2 ;il})D(x,y,'^)D(y 1 ,y 2 ;^) + 0(C 
Here the propagator D(x\, x;ip) is defined by 

dx£D(xi, x; if))W^ 2 \x, x 2 ) = 5(xi — x 2 ), 



while 



W^ m \x x ,...,x m ) 



5 m W 



8/j,(xi) . ..5/ji(x m )' 



(A9) 



(A10) 



(All) 



(cf. definition (|A1|) ). Restricting ourselves to the lowest order in the expansion flAjf ) (only 
including the saddle point terms proportional to () we obtain the form of the effective action 
used in the QSA approach (see, Eq.([l8|)). 
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FIGURES 



FIG. 1. Zero-temperature phase phase boundary in the Ej/Ec~vs-q x plane separating the 
Mott-insulating and superconducting states obtained from the quantum spherical approach (solid 
line), and the coarse-graining method (broken line). Note that there is no cusp at q x = since 
dEj(q x )/dq x \ qx= o = in the QSA. 

FIG. 2. Regular contribution of the frequency dependent conductivity at T = (real part) for 
several values of the dimensionless gap parameter 5 (that measures the distance from the critical 
point): a) 5 = 0.5, b) S = 0.1, c) 5 = 0.05 and d) 8 = 0.0. 

FIG. 3. Plot of the two parameter universal scaling function (real part) of the regular contri- 
bution to the frequency dependent conductivity in the quantum critical regime. 

FIG. 4. Same as in Fig. 3 for the imaginary part 
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Fig.l Quantum critical point 
T.K. Kopec & J.V. Jose 




Fig.2 Quantum critical point 
T.K. Kopec & J.V. Jose 




Fig.3 Quantum critical point 
T.K. Kopec & J.V. Jose 




Fig. 4 Quantum critical point 
T.K. Kopec & J.V. Jose 



